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Abstract 

In this paper, joint design of linear relay precoder and destination equalizer for dual-hop non- 
regenerative amplify-and-forward (AF) MIMO-OFDM systems under channel estimation errors is in- 
CZ3 , vestigated. Second order moments of channel estimation errors in the two hops are first deduced. Then 

based on the Bayesian framework, joint design of linear robust precoder at the relay and equalizer 
at the destination is proposed to minimize the total mean-square-error (MSE) of the output signal at 

t^^ I the destination. The optimal designs for both correlated and uncorrected channel estimation errors are 

!>■ 

If^ . considered. The relationship with existing algorithms is also disclosed. Simulation results show that the 



proposed robust designs outperform the design based on estimated channel state information only. 
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I. Introduction 

In order to enhance the coverage of base stations and quality of wireless links, dual-hop 
relaying is being considered to be one of the essential parts for future communication systems 
(e.g., LTE, IMT-Adanced, Winner Project). In dual-hop cooperative communication, relay nodes 
receive signal transmitted from a source and then forward it to the destination [[II, S. Roughly 
speaking, there are three different relay strategies: decode-and-forward (DF), compress-and- 
forward (CF) and amplify-and-forward (AF). Among them, AF strategy is the most preferable 
for practical systems due to its low complexity [[3][-[[7][. 

On the other hand, for wideband communication, multiple-input multiple-output (MIMO) 
orthogonal-frequency-division-multiplexing (OFDM) has gained a lot of attention in both in- 
dustrial and academic communities, due to its high spectral efficiency, spatial diversity and 
multiplexing gains [[8l-[[TT|. The combination of AF and MIMO-OFDM becomes an attractive 
option for enabling high-speed wireless multi-media services [[T2| . 

In the last decade, linear transceiver design for various systems has been extensively investi- 
gated because of its low implementation complexity and satisfactory performance [[HI, [[T3l . For 
linear transceiver design, minimum mean-square-error (MMSE) is one of the most important 
and frequently used criteria |[T4l - |[20l . For example, for point-to-point MIMO and MIMO- 
OFDM systems, linear MMSE transceiver design has been discussed in details in [[l4 | - [[T6l . 
Linear MMSE transceiver design for multiuser MIMO systems has been considered in [fTTI . For 
single carrier AF MIMO relay systems, linear MMSE precoder at the relay and equalizer at the 
destination are joint designed in [[l9l . Furthermore, the linear MMSE transceiver design for dual 
hop MIMO-OFDM relay systems is proposed in [[20l . 

In all the above works, channel state information (CSI) is assumed to be perfectly known. 
Unfortunately, in practical systems, CSI must be estimated and channel estimation errors are 
inevitable. When channel estimation errors exist, in general, two classes of robust designs can 
be employed: min-max and stochastic robust designs. If the distributions of channel estimation 
errors are known to be unbounded, stochastic robust design is preferred. Stochastic robust design 
includes probability-based design and Bayesian design. In this paper, we focus on Bayesian 
design, in which an averaged mean-square-error (MSE) performance is considered. Recently, 
Bayesian robust linear MMSE transceiver design under channel uncertainties has been addressed 
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for point-to-point MIMO systems ^^, [|23l and point-to-point MIMO-OFDM systems [|24ll . 

In this paper, we take a step further and consider the robust linear MMSE relay precoder and 
destination equalizer design for dual-hop AF MIMO-OFDM relay systems without considering 
direct link. For channel estimation in the two hops, both the linear minimum mean square error 
and maximum likelihood estimators are derived, based on which the second order moments 
of channel estimation errors are deduced. Using the Bayesian framework, channel estimation 
errors are taken into account in the transceiver design criterion. Then a general closed-form 
solution for the optimal transceiver is proposed. Both the uncorrected and correlated channel 
estimation errors are considered. The relationship between the proposed algorithm and several 
existing robust transceiver designs is revealed. Furthermore, simulation results demonstrate that 
the proposed robust algorithms provide an obvious advantage in terms of data mean-square-error 
(MSE) compared to the algorithm based on estimated CSI only. 

This paper is organized as follows. System model is presented in Section |nl Channel estimators 
and the corresponding covariance of channel estimation errors are derived in section Hill The 
optimization problem for transceiver design is formulated in Section |IVl In Section |Vl the general 
optimal closed-form solution for the transceiver design problem is proposed. Simulation results 
are given in Section |VI] and finally, conclusions are drawn in Section IVIIi 

The following notations are used throughout this paper. Boldface lowercase letters denote 
vectors, while boldface uppercase letters denote matrices. The notations Z^, Z^ and Z* denote 
the transpose, Hermitian and conjugate of the matrix Z, respectively, and Tr(Z) is the trace 
of the matrix Z. The symbol Im denotes the M x M identity matrix, while Oajxn denotes the 
M X A^ all zero matrix. The notation Z 2 is the Hermitian square root of the positive semi-definite 
matrix Z, such that Z = Z2Z2 and Z2 is a Hermitian matrix. The symbol E{.} represents the 
expectation operation. The operation vec(Z) stacks the columns of the matrix Z into a single 
vector. The symbol ® represents Kronecker product. The symbol a"*" means max{0,a}. The 
notation diag[A, B] denotes the block diagonal matrix with A and B as the diagonal elements. 

II. System Model 

In this paper, a dual-hop amplify-and-forward (AF) MIMO-OFDM relay cooperative com- 
munication system is considered, which consists of one source with Ng antennas, one relay 
with M/j receive antennas and Nj^ transmit antennas, and one destination with Md antennas, as 
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shown in Fig. \T\ At the first hop, the source transmits data to the relay, and the received signal 
Xfc at the relay on the k^^' subcarrier is 

Xfe = B.sr,kSk + ni,fc k = 0,1,-- -K -1, (1) 

where s^ is the data vector transmitted by the source with covariance matrix Rg^. = Efs^sf } 
on the k^^ subcarrier, and Rg^. can be an arbitrary covariance matrix. The matrix iisr,k is the 
MIMO channel between the source and relay on the k^^' subcarrier. The symbol rii ^. is the 
additive Gaussian noise with zero mean and covariance matrix Rnn. = (^ui^Mr ^^ '^he k^^ 
subcarrier. At the relay, for each subcarrrier, the received signal x^ is multiplied by a precoder 
matrix F^, under a power constraint ^^. Tr(FfcRxj.F^) < Pr where Rx^. = E{xfcX^} and Pr 
is the maximum transmit power. Then the resulting signal is transmitted to the destination. The 
received data y^ at the destination on the k'^^ subcarrier is 

Yk = Hj.rf,fcFfcHsr,fcS^. + H^rf^fcFfciiijfc + n2,fc, (2) 

where the symbol n2.A; is the additive Gaussian noise vector on the k^^ subcarrier at the second 
hop with zero mean and covariance matrix R^^ ^ = a'^J-M^- In order to guarantee the transmitted 
data Sfe can be recovered at the destination, it is assumed that Mji, Nji, and Md are greater than 
or equal to A^g [|6l. 

The signal x received at the relay and the signal y received at the destination in frequency 
domain can be compactly written as 

X = HsrS + ni, (3) 

y = UrdFUsrS + H,rfFni + na, (4) 

A r T T iT A r T T iT /c \ 

F^diag[Fo,---,F;,_i], (5b) 

Hsr = diag[lisrfl, iisr,i, ■ ■ ■ , iisr,K-i\, (5c) 

Urd = diag[Hrd,07 H^d,!, ■ ■ ■ , 'Hrd,K-l\, (5d) 

ni = Ko,n^i,--- ,n7,i^_i]^, (5e) 

^2 = [nJo, nji, ■ ■ ■ , nj;^_ J'^. (5f) 
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where 



Notice that in general the matrix F in (Hj) can be an arbitrary KNji x KMn matrix instead 
of a block diagonal matrix. This corresponds to mixing the data from different subcarriers 
at the relay, and is referred as subcarrier cooperative AF MIMO-OFDM systems EOl . It is 
obvious that when the number of subcarrier K is large, transceiver design for such systems 
needs very high complexity. On other hand, it has been shown in ll20ll that the low-complexity 
subcarrier independent AF MIMO-OFDM systems (i.e., the system considered in ([3]) and ^) 
only have a slight performance loss in terms of total data mean-square-error (MSE) compared 
to the subcarrier cooperative AF MIMO-OFDM systems. Therefore, in this paper, we focus on 
the more practical subcarrier independent AF MIMO-OFDM relay systems. 

III. Channel Estimation Error Modeling 

In practical systems, channel state information (CSI) is unknown and must be estimated. 
Here, we consider estimating the channels based on training sequence. Furthermore, the two 
frequency- selective MIMO channels between the source and relay, and that between the relay 
and destination are estimated independently. In our work, the source-relay channel is estimated at 
the relay, while the relay-destination channel is estimated at the destination. Then each channel 
estimation problem is a standard point-to-point MIMO-OFDM channel estimation. 

For point-to-point MIMO-OFDM systems, channels can be estimated in either frequency 
domain or time domain. The advantage of time domain over frequency domain channel estimation 
is that there are much fewer parameters to be estimated |[25l . Therefore, we focus on time domain 
channel estimation. Because the channels in the two hops are separately estimated in time domain, 
we will present the first hop channel estimation as an example and the same procedure can be 
applied to the second hop channel estimation. 

From the received signal model in frequency domain given by ([3]), the corresponding time 
domain signal is 

= (JP^ ® Im«)H,,(:F ® l^) (J^^ ® I^ )s + (:fH ® Ij,, )ni (6) 

^^ v '^ V ' ^^ V ■ ' 

where ^ is the discrete-Fourier-transform (DFT) matrix with dimension K x K. Based on the 
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properties of DFT matrix, it is proved in Appendix |A] that ^ can be rewritten as 



(DT®lM«)vec([?^ 



(0) 
sr 



n 



(Li-l) 



])+V, 



(7) 






where the matrices IV-^J are defined as 

K-l 



n^} = Y.^sr,ke^'^''\ £ = 0,1,---Li-1. 



(8) 



k=0 



It is obvious that H^^) is the t^^ tap of the multi-path MIMO channel between the source and 
relay in the time domain and Li is the length of the multi-path channel. The data matrix D is 
a block circular matrix as 

do di d/^„i 

d^--! do ■ • ■ • : d/^„2 



D 



d_ft:-Li+l ^K-Li+2 ^K-Li 



(9) 



where the element dj is expressed as 

I 

d 



^J^s^e^t^W, ^ = o,...,K-l. (10) 

Based on the signal model in ([7]), the linear minimum-mean-square-error (LMMSE) channel 
estimate is given by ll25ll 

with the corresponding MSE 

where Rchannci = 1E{^^^^3^^} is the prior information for channel covariance matrix. For 
uncorrelated channel taps, Rchannci = Kimnnd 'S)ImrNs and Achannci = diag[a ho, cr hi, ■■■ , o-/,^_J, 
where ahi is the variance of the /*^ channel tap ll24l . 

On the other hand, the channel in frequency domain and time domain has the following 



relationshi 



J 



Vec([H^,.,0 • • ■ iisr,K-l]) = ^{T L, ® ^MrNsMs 



(13) 



This relationship holds for both perfect CSI and estimated CSI. 
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where T j^^ is the first Li columns of J- . If the frequency domain channel estimate Hg^ ^ is 
computed according to (fT3l) . we have 



E{vec([AH,,,o ■■■ AH,,,x-i])vec^([AH,,,o ■■■ AH,,,,,_i])} 



(^Li ® ^MnNs) (Aeha 



channel ^ ^^S 



\f 



a-2(D*DT))~l ®\mA^L, ® ^MnNsfK. (14) 






where AH^j. ^ = H^j. ^ — iisr k- 

In case there is no prior information on Rchannch we can assign uninformative prior to ^^r, that 
is, cT/io : CT/ii ; ■ ■ ■ ) CT/iL-i approach infinity ll26l . In this case, Rj^hinnci ~^ 0' ^^^^ then the channel 
estimator (fTTI) and estimation MSE (fT2l) reduce to that of maximum likelihood (ML) estimation 
nisi P. 179]. 

Taking the MrNs x MrNs block diagonal elements from ((14)) gives 

/Li-lLi-l 

E{vec(AH,,,fc)vec«(AH,,,,.)} = E E (^"'^'^''"''^K< 

\^2=0 £i=0 

where $£[ ^^ i^ the A^5 x A^^ matrix taken from the following partition of ^"^ 



/2' 



I 



Mr- 



(15) 



$* 



^0,0 



^ai 



^0,Li-l 



$ 



Li-1,0 






(16) 



Furthermore, based on (fTSi) . for an arbitrary square matrix R, it is proved in Appendix |B] that 

E{AH..,,RA<J = Tr [rJ^ Y. {^~''^'^''-''\^Z,,f)] Im«- (17) 



=0 fi=0 



A similar result holds for the second hop. In particular, denoting the relationship between the 
true value and estimate of the second hop channel as 



Hr-Wfc — Hrd.fc + AHrdj 



A; = 0,--- ,/s:-i, 



we have the following property 

E{AH.,,,RAH^,,J = Tr ( Rj^ 5^ (e-^1r'=(^i-^^)($- 



L2-1L2-1 



rd \T 

12) 



lAfr 



(18) 



(19) 



=0 £2=0 



where L2 is the length of the second hop channel in time domain. Furthermore, as the two 
channels are estimated independently, AHsf,,fc and AHf^i ^ are independent. 
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IV. Transceiver Design Problem Formulation 

At the destination, a linear equalizer G^ is adopted for each subcarrier to detect the transmitted 
data Sk (see Fig. [U). The problem is how to design the linear precoder matrix F^ at the relay 
and the linear equalizer G^ at the destination to minimize the MSE of the received data at the 
destination: 

MSEfc(Ffc, Gfc) = E{Tr ((G^yfe - Sfc)(Gfcy, - s^)")}, (20) 

where the expectation is taken with respect to s^, AH^j. fc, AH^rf fc, rii ^ and n2^k- Since s^, rii ^ 
and n2,fe are independent, the MSE expression (|20l) can be written as 

MSEfc(Ffc,GJ 
= E{||(GfcHrrf^fcFfcHsr,fc — lAfs)sfc + GfcHj.rf^fcFfcni^fc + Gfcn2,fc|| } 

= lE,AHsr,fc,AHrd,fe{Tr((GfcHj.rf^fcFfcHsr,fc — iNs)'^Sk{Gkiird,kFkiisr,k — ^Ns) )} 

+ Eah.,, jTr {{GkH.ra,kFk)Rn,AGkiird,kFkf)} + Tr(GfcR„, ,G^) 

= ^AHsr,kA'H-rd^k{'^^ {i^kiird,kFkiisr,k)^Sk{Gkiird,kFkiisr,k) )} 

+ Tr (GfcEAH,d,fc{Hrd,feFfcR'ni,fcFfH^^_fc}Gf ) 

— Tr \Ils,k{G kiird,kF kiisr,k) j — Tr {Gkiird,kFkiisr,k^s,kj 

+ Tr(R,J + Tr(G,R„,,,G«). (21) 

Because AH^r,/; and AH^.^ ^ are independent, the first term of MSE^ is 

^AH,r,kAiird,k{'^''^ {{Gkiird,kFkiisr,k)'Rsk{Gkiird,kFkiisr,k) )} 

= Tr (GfcEAH,d,fe {Hrd,fcFfcEAH,,,fc{Hsr.,fcRsfcH^r,fc}FfeH^d,fc} G^) . 
For the inner expectation, the following equation holds 

lEAH^^,fc{Hsr,fcRsfeH^^^;,} = EAHs,,,fe{(Hsr-,/c + AHsr,fc)Rsfe (Hsr-,fc + AHgr^k) } 

= TlCRsk^ sr,k)'^MR + ^sr,k^Skiisr,k = FEfc, 

where based on ([17]) the matrix ^sr,k is defined as 

^sr,k = E E (e--'^^^^-^^) in.f) ■ (24) 



^1=0 ^2=0 



(22) 



(23) 
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Applying (|23l) and the corresponding result for AUrd^k to (l22l) . the first term of MSE^ becomes 

Tr (GfeEAH,d,fe {H.rd,fcFfcEAH,,,fe{H^r,fcRsfeH^r,fc}FfeH^d,fc} G^) 
= TT{Gk{TT{FkUkFf^rd,k)lMn + H.rf,fcFfcnfcF^H,%)G») , (25) 

where the matrix ^rd,k is defined as 

^r.,k = y: y: {e-^'^^'^'^-^^KKif) ■ (26) 

£1=0 £2=0 

Similarly, the second term of MSE^ in (|2T1) can be simplified as 

Tr (GfcEAH^j,fc{Hrd,fcFfcR„j_feFfcH^^;,}G^,) 

= Tr(Gfc \Tr{FkIlni^k^k ^rd,k)^Mn + Hrd,fcFfcR„j_fcF^ H^^^^ j Gf, )■ (27) 

Based on ([25]) and ([27]), the MSE^ (ED) equals to 

MSEfc(Ffc, Gfc) = Tr (^Gfc(Hrd,fcFfcRxfcFf HJ:'^ ;, + Kfc)Gf j - Tr (^Rg^H^^ ;,F^H^^ ^.G^ ^ 

— Tr ( GfcHrd,fcFfcHsr,A;Rsfc j + Tr(Rsj^) (28) 

where 

Rx, =n, + <lM, (29) 

Kfc = (Tr(F,Rx,F«*,,,fc) + Ol,,, 

^Vk^MD- (30) 

Notice that the matrix Rx^ is the correlation matrix of the receive signal x^ on the k^^ subcarrier 
at the relay. 

Subject to the transmit power constraint at the relay, the joint design of precoder at the relay 
and equalizer at the destination that minimizes the total MSE of the output data at the destination 
can be formulated as the following optimization problem 

min VMSEfc(Ffc,Gfc) 

k 

s.t. 5^Tr(FfcRx,Ff)<P,. (31) 

k 

Remark 1: In this paper, the relay estimates the source-relay channel and the destination 
estimates the relay-destination channel. The precoder F^ and G^ can be designed at the relay 
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or at the destination. Sharing channel estimation between the relay and the destination is un- 
avoidable. However, when channel is varying slowly, and the channel estimation feedback occurs 
infrequently, the errors in feedback can be negligible. 

V. Proposed Closed-Form Solution 

In this section, we will derive a closed-form solution for the optimization problem (|3T| ). In 
order to facilitate the analysis, the optimization problem (|3T| ) is rewritten as 

min VMSEfc(Ffc,G,) 

Fk,Gk,Pr,k ""^ 
fe 

s.t. Tr(FfcRx,Ff ) < Pr,k, k = 0,--- ,K-1 

Y,Pr,k<Pr. (32) 

k 

with the physical meaning of Pr^k being the maximum allocated power over the k^^ subcarrier. 
The Lagrangian function of the optimization problem (l32l) is 

£(F,, Gfc, Pr,k) = Y, MSEfc(F,, Gfc) + J^ 7fc(Tr(FfcR,,Ff ) - P,,,.) + 7(5^ P,,,. - P,) 

(33) 
where the positive scalars 7^ and 7 are the Lagrange multipliers. Differentiating (l33l) with respect 
to Fjfc, Gfc and P^ fc, and setting the corresponding results to zero, the Karush-Kuhn-Tucker (KKT) 
conditions of the optimization problem (|32l ) are given by ||27l 

Gfc(Hrd,fcFfcRxfcF^H^^;, + Kfc) = Rsj^(Hrd,A;FA;Hsr,fc) , (34a) 

Hr.d,fcGfcGfcHrd,A:FA;Rxfc + (Tr(GfcGfc)*r.d,fc + 7fc)F/cRx^ = (Hsj,,fcRsfeGfcHrd,fcj , (34b) 

7fc(Tr(FfcRx,F^) - P,,^) = 0, (34c) 

7fc>0, k = 0,---,K-l, (34d) 

7(^P,,fc-P,) = 0, (34e) 

k 

7o = 7i = ■ ■ ■ = iK-i = 7, (34f) 

Tr(FfcRx,Ff ) < P,,fc, (34g) 

5^Pr,fc<Pr. (34h) 

k 
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It is obvious that the objective function and constraints of (|32l ) are continuously differentiable. 
Furthermore, it is easy to see that solutions of the optimization problem (|32|) satisfy the regularity 
condition, i.e., Abadie constraint qualification (ACQ), because linear independence constraint 
qualification (LICQ) can be proved Il28l . Based on these facts, the KKT conditions are the 
necessary conditionso From KKT conditions, we can have the following two useful properties 
which can help us to find the optimal solution. 
Property 1: It is proved in Appendix O that for any F^ satisfying the KKT conditions (|34a|) - 



(|34el) . the power constraints ( |34g[ ) and (I34h|) must occur on the boundaries 

Tr(FfcR,,F^) = Pr,k, (35) 

Y,Pr,k = Pr- (36) 

k 

Furthermore, the corresponding G^ satisfies 

Tr(GfcG») = 7fcP,,fcK^. (37) 

Property 2: Define the matrices Ut^., Vtj. , At^., U©^, and A©^. based on singular value 
decomposition (SVD) as 

iPr,k^rd,k + al,lNn)~^B-fd,k '^rd,kiPr,k^rd , k + ^IjV^)'^^ = Uq^Aq^UI^, (38) 

=@k 

R:|Hs.,fcRs, =Ut,At,V?^, (39) 

with elements of the diagonal matrix Atj. and A^^ arranged in decreasing order. Then with 
KKT conditions (I34al) and (I34bl) . it is proved in Appendix |D] that the optimal precoder F^ and 
equalizer G^ for the optimization problem (l32l) are in the forms of 



Fk = {Pr,k^rd,k + <I;vJ-^Ue„,,AF,U?^,,p^Rx,^ (40) 

G, = Yrr,,,,AG,Vl^^jPr,,^rd,k + <I^J-iH«,„ (41) 

^Notice that the solution Fo = ■ ■ • = Fk-i ~ and Go — ■ ■ ■ = Gk-i ~ also satisfies the KKT conditions, but this 
solution is meaningless as no signal can be transmitted 1141 . 
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where A-p^^ and Aq^ are to be determined. The matrix Ut^^.p^ and 'VTk,Pk ^^ the first pk columns 
of Uxfe and Vtj^, respectively, and pk = Rank(ATj^). Similarly, U^j^^q^ is the first g^ columns 
of Ue^^, and q^ = Rank(A@^). 



Right multiplying both sides of (I34al) with G^ and left multiplying both sides of (I34bl) with 
F^, and making use of (|40l) and (BTl) . the first two KKT conditions become 



AG,Ae,AF,A«^A0,A«^ 



%^G,Ae,A« 



A«, Ap,. + ^AP Ap 



-^F^^k^Gk^Gk^'^&k^Fk ^ -^^Ffc-^Ffc 



(AG,Ae,AF,ATj^, 



(At.Ag.Aq.AfJ 



(42) 
(43) 



where the matrix A©^ is the qk x q^. principal submatrix of A©^ . Similarly, At^ is the pk x p^ 
principal submatrix of Atj.. In this paper, we consider AF MIMO-OFDM relay systems, the 
matrices A^^ and Aq^ can be of arbitrary dimension instead of the square matrices considered 
in [[T4| and [[22l|. It should be noticed that as the optimization problem (l32l) is not a convex 
problem, the KKT conditions are only necessary conditions. That is, there are many solutions 
that will satisfy the KKT conditions. To identify the optimal solution, we need an additional 
information which is presented in the following Property 3. 

Property 3: Putting the results of Property 1 and Property 2 into the optimization problem 
(|32|) . based on majorization theory, it is proved in Appendix |E] that the optimal A^^ and Aq^ 
have the following diagonal structure 



(44) 



^Ffe,opt — 


Ap^^opt 

Ogk-Nk,Nk 


^qk-Nk,Pk~Nk 


Acfejopt = 


Acfc.opt 


ONk,qk-Nk 




_ Op, -TV,, TV, 


%k-Nk,qk~Nk 



(45) 



where AF^.,opt and AG^.,opt are two A^^ x Nk diagonal matrices to be determined, and Nk = 
mm{pk,qk)- 

Substituting (|44l) and (l45l) into (|42l) and (l43l) . and noticing that all matrices are diagonal. 
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^Ffc,opt and Acfc.opt can be easily solved to be 



A 



Ffc.opt 



A, 



Gk,opt 







^fcAg, 



7fc 



VkC^l, 



A0^,At. - 4^ A 



2 ^^©fc 



0", 



n2 



Ae^^ 



(46) 



(47) 



where the matrices A^^. and A©^ are the principal sub-matrices of At^. and Aq^. with dimension 
Nk X Nk, and Nt = min{rank(A©J,rank(ATj}. The matrices U^^^n^^ ^Tk,Nk and Uqi^^n^ 
are the first N^ columns of Ut^, Vt^^ and \J@^, respectively. 

In the general solution P6l)-(|47]). Pr^k, Vk and 7^ are unknown. However notice that from ( 1351 ) 
and (l37l) in Property 1 , the optimal precoder and equalizer should simultaneously satisfy 

Tr(Ffc^optR-XfcFfcQpJ = P,.fe, (48) 

Tr(Gfc,optG«„pJ = 7/cPr,fc/<. (49) 

Substituting (|44|)-(|471) into (l48l) and ( |49l) . it can be straightforwardly shown that rjk and 7^ can 
be expressed as functions of P^.fc 

b3,kPr,k 



Vk 



Ik 



Pr,kbl,k + &l,fc&4,fc — &2,/t&3,fc 
b-i,k(^n2iP'r,kbl,k + ^l,fe^4,fc " ^2,fc^3,fc) 



(-fr,fc + &4,fc)^-fr,fc 

where 61 fc, 62,^, ^3,^ and 64 ^ are defined as 

6i,fc ^ Tr(U|^,^^(P,,fc*,,,fc + <I^J"iU©„^,AT,AejAi,,), 

&2,fc = Tr(U0^_^^(Pr,fc*rd,fc + O-naW)" Ue^^^Ar^^ Aq^ Aj^fc), 

63,^ ^ Tr(AT,AejAi,,), 
&4,fc = Tr(A@ Ai^fc), 



(50) 
(51) 

(52a) 
(52b) 
(52c) 
(52d) 



and Ai fc is a diagonal selection matrix with diagonal elements being 1 or 0, and serves to replace 
the operation '+'. Combining all the results in this section, we have the following summary. 
Summary: The optimal precoder F/,. opt and equalizer GA,-,opt are 
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A;,opt 



Gfc,opt = VTfc.ATj^AGfe.optUg,^ jVfc(-Pr,fc*rd,fc + (^ni^NR) ^H^d,fc, 



(53) 
(54) 



where 



^Ffc,opt 



^Gk,opt 




7fc 



^' -AZAn 



%f^: 



"2 



^n2 



A ^ 



(55) 



(56) 



with r]k and 7^ given by (|50l)-(l52]). 

From the above summary, it is obvious that the problem of finding optimal precoder and 
equalizer reduces to computing Pr^k, and it can be solved based on (IST]) and the following two 
constraints (i.e., (l34fl) and (l36l) ) 



7o 



r,fc 



7i^-i, 
P,. 



(57) 
(58) 



In the following subsections, we will discuss how to compute Pr^k- 

Remark 2: When both channels in the two hops are flat-fading channels, the considered 
system reduces to single-carrier AF MIMO relay system. It should be noticed that for single- 
carrier MIMO relay systems, there is no need to consider power allocation among subcarriers, 
and we can set Pr,fe = Pr- In this case, the proposed closed-form solution is exactly the optimal 
solution for the robust transceiver design in flat-fading channel. Furthermore, when the CSI in 
the two hops are perfectly known, the derived solution reduces to the optimal solution proposed 
in|fT9l. 

Remark 3: Notice that when the source-relay link is noiseless and the first hop channel 
is an identity matrix, the closed-form solution can be simplified to the optimal robust linear 
MMSE transceiver for point-to-point MIMO-OFDM systems Il24ll . Moreover, if single carrier 
transmission is employed, the closed-form solution further reduces to the optimal point-to-point 
MIMO robust LMMSE transceiver It22l . 
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Remark 4: The complexity of the proposed algorithm is dominated matrix decomposition, 
matrix multiplication and matrix inversion. The complexity of each of these matrix operations is 
known to be 0{n^), where n is the matrix dimension Il29ll . So the complexity of our proposed 
algorithm is 0{m^), where m = maxjAf^), Nji, Nd, Ns}. 

A. Uncorrelated Channel Estimation Error 

When the channel estimation errors are uncorrelated (for example, by using training sequences 
that are white in both time and space dimensions), the following condition must be satisfied filOil . 
Il30l-|l32l 

DD^ocI^.i,. (59) 

Then according to (fT4l) . we have ^ sr,k = Xl^ ^T i 1^ ^ '^Ns- Similarly, for the second hop, 
we also have 

^rd,k OC Iat^ = 5rd,k^Na-, (60) 

where the specific form of 5rd,k can be easily derived based on (|26|) . 
Putting (l60l) into the left hand side of (l38l) . the expression becomes 

{Pr,k^rd,k + (jlj.Nn)~^'tl'^d,k^rd,k{Pr,k^rd,k + <I7V«)~^ 

Hrd.feHrdt. (61) 



Pr,k^rd,k + '7^2 



Applying eigen-decomposition H^^^Hrd.fc = UHfcAHj,UH and comparing with the right hand 
side of (|38l) . we have 

U0, = Uh,, Ae, = -p^- --^Ah,. (62) 

[-t^r,kOrd,k + <^n2) 



Substituting (1621) into (1511) . 7^ reduces to 



< (Tr(AT,AHjAi,,)) 

Ik = ^— ^ ^, (63) 

(P,,fc (1 + 6rd,kTriA^[Ai^k)) + alTi{A^[AiS 

where Ah^, is the N^ x N^ principal submatrix of Ah^, • 
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With (|63l) and the facts that ^^ Pr^k = Pr and 7o = ■ ■ ■ = 'Jk-i, Pr,k can be straightforwardly 
computed to be 



1 



r -1 



(7^ Tr(AT,A„'Aifc) a2Tr(A„Ai^ 

7 1 + 6rd,kM^uA^,k) l + VfcTr(AH,Ai,fc) 
where 7 equals to 

-,^<|^ TV(A..A„!a„) V/f,,^^ <TV(A-A„) y ^^^^ 

fc l + (5,.rf,fcTr(AH^Ai,fc.)y / V fc l + ^rd,fcTr(AH^Ai,fc)y 

5. Correlated Channel Estimation Error 

Due to limited length of training sequence, DD^ oc I may not be possible to achieve [f30l . 
In this case, the channel estimation errors are correlated, and ^rd,k 9^ I- From (|38l) . it can be 
seen that the relationship between A©j, and P^.fc cannot be expressed in a closed-form . Then the 
solution for P^.fc cannot be directly obtained. However, notice that when Pr, k^rain{^ rd,k) ^ cr^2' 
where Anun(Z) denotes the minimum eigenvalue of Z, we have 

Pr,k^rd,k + (^ui^Nfi ~ Pr,k^rd,k- (66) 

This situation occurs at high SNR in the second hop, and we term this high SNR approximation 
(HSA). On the other hand, when Pr,k^mm{^ rd,k) ^ o"^2 cannot be guaranteed, we have 

Pr,k^rd,k + O'nalA'fl ~ (-Pr,fcAmax(*rd,fc) + O^nJlAffl, (67) 

and it is termed spectral approximation (SPA). For spectral approximation, ^rd is replaced by 
Amax(^rd)I, and from the MSB formulation in (l28l) . it is obvious that the resultant expression 
forms an upper-bound to the original MSB. Notice that when the training sequences are close 
to white sequence [[34l . Il35l . the eigenvalue spread of ^rd is small, and SPA is a good approx- 
imation. In the following, computations of Pf.,fc under different approximations are detailed. 

1) High SNR Approximation (HSA) at the second hop: Based on HSA, the left hand side of 
(l38l) becomes 

{Pr,k'^rd,k + al^lN^y^B.fa,k^rd,k{Pr,k^rd,k + <l7Vfi)"^ 

1 _H . _i 

2 "LjH "LJ ,Tr 2 



P 



VTV 2 TLTn TT Vlf 2 

^rd,k"-rd,k "-rd,k^rd,k 

Ur,Ar,U», (68) 



r,k 
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where the second equality is based on SVD. Comparing the right hand side of (|681) with that of 
(l38l) . we directly have 

Ue, =Ur„ Ae, = -^Ar.. (69) 

Substituting (|69l ) into (|5T1) . 7^ reduces to a simpler form 

1 C3,A;Cr^2(ci,fc + Ci,fcC4,fc -C2,fcC3,fc) 
^ . ' 

A 

=Xk 

where ci^k, C2,k, c^^k and 04,^ are defined as 

ci,fc ^ Tr(U«^,^^*;iUr„^,AT,Arf Ai,fc), (71a) 

C2,fc ^ Tr(U^^,^^*;,^,Ur„^,Ar;Ai,fc), (71b) 

C3,fc = Tr(AT,ArfAi,fc), (71c) 

C4,fc = Tr(Ar'Ai,fc), (7 Id) 

with the diagonal matrix Ap^ being the A^^ x Nk principal submatrix of Ar^., and the matrix 
Urfe.ATfe consists of the first A^^ columns of the matrix Ur^. Together with the facts that ^^ P^.fc = 
Pr and 7o = ■ ■ ■ = 7^--!, Pr,k can be solved as 

^ v^ (72) 



Y.k\^ 

2) Spectral Approximation (SPA): With SPA, the left hand side of (|38l) becomes 



{Pr,k^rd,k + Crri2-'-AfH) ^ ^rd,k^rd,k{Pr,k^ rd,k + <^n2-'-Affl) ^ 



-* r,fe^maxl^^rd,/cj T 0", 



2 

n2 



Hr-d.fcHr-WA;- (73) 



Comparing (1731) to (|6T1) . it is obvious that the problem becomes exactly the same as that discussed 
for uncorrelated channel estimation errors. Therefore, the allocated power to the k^^ subcarrier 
Pr,k can be calculated by (l64l) but with 5rd,k replaced by Xma.xi^rd,k)- 

Notice that the selection criterion Pr,kXmmi^rd,k) ^ 0-^2 involves the parameter of interest 
Pr^k- In practice, we can replace Pr,fc in the criterion by Pr/ K. Then the proposed algorithm 
can be summarized as using HSA when (-Pr/A^)Amin(*rd,fc)/o"^2 — ^' otherwise we use SPA, 
where T is a threshold. 
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Remark 5: Our design is valid for both cases whether the source has channel state information 
(CSI) or not, since our design is suitable for any correlation matrix Rg ^ which is determined by 
the precoder at the source. The precoder at the source can be easily designed based on different 
criteria such as zero forcing (ZF), capacity maximization (CM) and minimum-mean-square-error 
(MMSE) using the first hop channel information. Once the precoder at the source is fixed, the 
optimal forward matrix at the relay and the equalizer at destination can be designed using the 
proposed design. 

VI. Simulation Results and Discussions 

In this section, we investigate the performance of the proposed algorithms. For the purpose of 
comparison, the algorithm based on estimated channel only (without taking the channel errors 
into account) is also simulated. An AF MIMO-OFDM relay system where the source, relay 
and destination are equipped with same number of antennas, Ns = Mr = Nr = Mr, = 2 
is considered. The number of subcarriers K is set to be 64, and the length of the multi-path 
channels in both hops is 5, and is denoted as L. The channel response is generated according 
to the HIPERLAN/2 standard ^^. The signal-to-noise ratio (SNR) of the first hop is defined 
as Es/Ni = Tr(Rs)/(-ft'Mijcr^J, and is fixed as 30dB. At the source, on each subcarrier, two 
independent data streams are transmitted by two antennas at the same power, and QPSK is used 
as the modulation scheme. The SNR at the second hop is defined as Er/N2 = Pr/{KMcia'^J. 
In the figures, MSB is referred to total simulated MSB over all subcarriers normalized by K. 

Based on the definition of D in Q, DD^ is a block circular matrix. In the following, only 
the effect of spatial correlation of training sequence is demonstrated, and the training is white 
in time dimension. In this case, DD^ is a block diagonal matrix, and can be written as DD^ = 
II ® J2i djdp, where J2i djdf /fC is the spatial correlation matrix of the training sequence. 
Furthermore, the widely used exponential correlation model is adopted to denote the spatial 
correlation matrix [[22l|. Il23l . and therefore we have 

1 a 
a 1 

It is assumed that the same training sequence is used for channel estimation in the two hops. 



DD^ = lr^K 



(74) 
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Based on the definition of ^sr,k and '^rd,k in (l24l) and ( |26l ). and together with (1741 ), we have 

-1 



_ L 2 

^sr,fc ^rdjfc "E^^e 

A 



1 a 
a 1 



(75) 



where a^ is the noise variance during channel estimation. 

First, we investigate the performance of the proposed robust algorithm when channel estimation 
errors are uncorrected, which corresponds to a = in (1751) . Fig. [21 shows the MSE of the received 
signal at the destination with different a^. It can be seen that the performance of the proposed 
robust algorithm is always better than that of the algorithm based on estimated CSI only, as long 
as a"^ is not zero. Furthermore, the performance improvement of the proposed robust algorithm 
over the algorithm based on only estimated CSI enlarges when a'^ increases. Fig. [3] shows 
the corresponding performance under correlated channel estimation errors (a = 0.4), and the 
threshold T is set to be 10. It can be seen that a similar conclusion can be drawn as in Fig. [2l 

Fig. |4] shows the MSE of the received signal at the destination for HSA, SPA and the proposed 
algorithm that switches between the two, when a = 0.4. It is clear that HSA performs better than 
SPA at high SNR region, as at high SNR, (l66l) in HSA becomes equality. On the other hand, 
SPA performs better than HSA at low SNR region. The proposed robust algorithm combines the 
benefits of both HSA and SPA. 

Fig. [5] shows the MSE of the output data at the destination for both proposed robust algorithm 
and the algorithm based on estimated CSI only, with different a and Er/N2 = 25dB . It can be 
seen that although performance degradation is observed for both algorithms when a increases, 
the proposed robust algorithm shows a significant improvement over the algorithm based on 
estimated CSI only. Furthermore, as a = gives the best data MSE performance, it demonstrates 
that white sequence is prefered in channel estimation. 

Finally, Fig. [6] shows the bit error rates (BER) of the output data at the destination for different 
a^, when a = 0.5. It can be seen that the BER performance is consistent with MSE performance 
in Fig. [2] and Fig. [3l 

VII. Conclusions 

In this paper, linear robust relay precoder and destination equalizer were jointly designed for 
AF MIMO-OFDM relay systems based on MMSE criterion. The linear channel estimators and the 
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corresponding MSE expressions were first derived. Then a general solution for optimal precoder 
and equalizer was proposed. When the channel estimation errors are uncorrected, the optimal 
solution is in closed-form, and it includes several existing transceiver design results as special 
cases. On the other hand, when channel estimation errors are correlated, a practical algorithm was 
introduced. Simulation results showed that the proposed algorithms offer significant performance 
improvements over the algorithm based on estimated CSI only. 



Appendix A 
Proof of & 

Based on the characteristics of DFT operation, the matrix l-Lsr defined in ([61) is a KMr x KNs 
block circulant matrix 

' u^^) ■■■ nf,}-''^ ni^.^-'^ ■■■ ?^«" 
?^« ?^1? ■■■ 7^1^-1) ■■■ 7^1?) 



rLf, 



"sr ' "ST 



(76) 



■■•0 T^i^"^) uf,}^'^^ nf,'-^'^ ■■■ nf) 

where the element 'H}{') is defined in ([8]). It is obvious that l-O-^}. is the t^ tap of the multi-path 
MIMO channels between the source and relay in the time domain and Li is the length of the 
multi-path channel. 

On the other hand, based on the definition of d in ^, we have 



do 



di 



dif-i 



From (1761) and (1771) . by straightforward computation, the signal model given in Q can be 
reformulated as 



r = -Hsrd + V = Yec{[Hf^ ■ ■ ■ Ui";'^'^]!)) + v 
where the matrix D is defined in Q. 



(78) 
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Appendix B 
Proof of (fTTI) 

For the expectation of the following product 

S = E{QRW"} (79) 

where Q and W are two M x N random matrices with compatible dimension to R, the (i, j)*^ 
element of S is 

E(^, j) = E{Q(z, :)RW(j, :)«} = J2J2 ^^Q^^' ^)^(^' ^) W(^' ^)*}- ^^O) 

t k 

If the two random matrices Q and W satisfy 

E{vec(Q)vec^(W)} = A ® B, (81) 

where A is a NxN matrix while B is a MxM matrix, then we have the equality E{Q(ii, ji)W(z2, J2)*} 
B(zi, Z2)A(ji, J2)- As Q(i,t) and W(j, A;) are scalars, (l80l) can be further written as 

Finally, writing (f82l) back to matrix form, we have [[36ll 

S = BTr(RA'^). (83) 

Notice that this conclusion is independent of the matrix variate distributions of Q and W, but 
only determined by their second order moments. Putting A = J2e ^=0 12e =0 (e~-'"^'^'^^^~^^^$^[ gj, 
B = Ia/^ and Q = W = AHsr^k , into ([83]), we have (flTl) . 

Appendix C 
Proof of Property 1 

Right multiplying both sides of (I34al) with Gf, the following equality holds 

Gfc(Hr(i,fcFfcRxj.F^H^^^^, + Kfc)G^ = Rs;.(Hrd,fcFfcHsr,fc) G^. (84) 

Left multiplying (I34bl) with F^, we have 

Ffc H^d.feG;, GfcHrd,fcFfcRxj, + F^, Tr(GfcGfc )*r'd,fcFfcRxfe + 7fcF^ F^Rx^^ 
— ^k ( Hsr.fcR-sfeGfcHrd^fc j . (85) 
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After taking the traces of both sides of ( f84l ) and ([85] ) and with the fact that the traces of their 
righthand sides are equivalent, i.e., TT{Rs^{'Hrd,kFkilsr,k)^Gf) = TT{F^{H,r,k'RskGktlrd,k)^), 
we directly have 

Tr(Gfc(Hrrf,fcFfcRxfcFfc li^d,k + ^k)Gk ) 

= Tr(F«H,",^,G^G,H,,,,F,R.J + 7fcTr(Ff F^R^J + Tr(G,Gf)Tr(Ff *,,,,F,R.J. 

(86) 

By the property of trace operator, Tr(Gfc(H,,,fcFfcR,,FHHH ,)Gf) = Tr(FHHH _,G^GfcH,,,fcFfcR,J, 
and (|86l) reduces to 

Tr(GfcKfcG^) = Tr(G,G^)Tr(F^*,rf,fcFfcR,J + 7feTr(F«F,R,J. (87) 

On the other hand, based on the definition of K^ in (l30l) . Tr(GfcK/cG^) can be also expressed 
as 

Tr(GfcKfcGf ) = Tr(GfcGf )Tr(F^*,rf,fcFfcRxJ + Tr(GfcR„,,fcG«). (88) 

Comparing (l87l) with (f88l) . it can be concluded that 

Tr(GfcR„,,,,Gf ) = 7feTr(FfcR,,Ff ). (89) 

Putting (EH) into dMc]), we have Tr(GfcR„2,fcG^) - -fkPr,k = 0. As R^^ , = a^Mo^ it is 
straightforward that 

<Tr(GfcG^) = 'ykPr,k. (90) 

Furthermore, based on the fact 70 = 7i = • • • = 7k- 1 = 7 and taking summation of both 
sides of (l90l) . the following equation holds 



^<Tr(G,Gf) = 7E^'-.^- (91) 

k k 

Putting dn]) into (I34el) . we have 

5^a^^Tr(G,G^)-7P. = 0, (92) 

k 

and it follows that 

7. = 7 = <^^^^^P^. (93) 



Pr 



June 1, 2010 DRAFT 



23 

Since for the optimal equalizer G^, ^^ Tr(Gfc^optG^QpJ 7^ 0, it can be concluded that 7^ 7^ 0. 
In order to have (|34c|) satisfied, we must have 

IV(Ffc,,ptRx,FH„pJ = Pr,k. (94) 

Furthermore, as 7 7^ 0, based on (|34el ). it is also concluded that 

Y,Pr,k = Pr- (95) 

Finally, (l90l) constitutes the second part of the Property 1. 

Appendix D 
Proof of Property 2 

Defining a full rank Hermitian matrix M^ = Pr,k^rd,k + o'n2^Nii, then for an arbitrary Nr x Nr 
matrix F^, it can be written as 

1 ,,1 



F,. = M,^Ue,SF,U5^Rx,^ (96) 

1 1 
where the inner matrix Sp^ equals to Sp^^ = U@ M^FfcRxfeUT^- 

Putting (|96| ) into (I34al) . and with the following definitions (the same as the definitions in (1381) 

anddH) 

M,:^H»,,H,,,,,M7 = Ue,A0,U|^., (97) 

_i ^ 

■Rx.fcHsr.fcR-s,*: = Ux^ATfeVrj,^, (98) 

the equalizer G/,. can be reformulated as 

Gfc = HskC^rd^kF kiisr,k) (Hrd.fcFfcRxfcF^ H^^ ;, + rjklMj^) 

= (Rxfe" Hsr^fcRs J (R:J^F^H^^;,H,,rf^fcFfcR:^,^ + ?7fclA/^)" K^^Fi^H^^j^ 

= Vt, A?^(SP^ Ae,EF, + r^klM^y'll^^^ U|^M,-^H« ,„ (99) 



where the second equality is due to the matrix inversion lemma. 

_i 
Putting (l90l) from Appendix O into (I34bl) . after multiplying both sides of (I34bl) with M^. % 

we have 

M^^Hj^^,G^GfcH,d,fcFfcRi + MfF^Ri^ = M^ (llsr,k'Rs,kGkilrd,kY Ti'l (100) 
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Then substituting F^. in ([961) and G^ in (|99l) into (IIOOI ). we have 

Sf = (Ae,S»^So,Ae, + ^I^J"^(At,Sg,A0J^. (101) 

Since At^. and Aq^. are rectangular diagonal matrices (denoting their ranks by pk and qk 
respectively), based on (IIOII) . it can be concluded that Sp^. has the following form 



^Ffc 







(102) 



J NrXMr 

where Ap^ is of dimension qk x pk and to be determined. Furthermore, putting (11021) into the 
definition of Sq^^ in (|99| ). we have 



Sg. 



^G. 





(103) 



J NsxMd 

where Aq^ is of dimension pk x qk, and to be determined. Substituting (11021 ) and (11031 ) into 
(|96l ) and (199]), it can be concluded that 



where 






Ag, = A;^(AP^A0,Af, +%IpJ U^^ 



(104) 
(105) 

(106) 



and Atj. is the p^ x pk principal submatrix of At^.. 

Appendix E 
Proof of Property 3 

Taking the trace of both sides of (|42l) and (l43l) . and noticing that the resultant two equations 
are the same, it is obvious that 



Tr(AG,Ae,A« 



-Tr(A«Ap,). 



'' %< ^ '' " 



On the other hand, substituting (11051) into (l90l) in Appendix O we have 



Tr(AG,Ae,A^J 



7fc_ 

2 



0", 



-'r,fc- 



(107) 



(108) 



"2 
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Comparing (11071) and (I108D . it follows that 

-Tr(A»^ AfJ = Pr,k. (109) 

For the objective function in the optimization problem (|32|) . substituting (l40l) and (|4T|) into 
the MSB expression in (l28l) . the MSB on the k^'^ subcarrier can be written as 

MSEfc(F,,Gfc) 
= Tr(Ai^(lAi^Ae,AF, + I,J-^) + Tr(R,,) - Tr(R,,H«,,R;^iH,,,,R,J, (110) 

Vk ^ ' 

A 
=Cfe 

where c^ is a constant part independent of F^. Therefore, based on (11091) and (IllOL the 
optimization problem (|32l) becomes as 

min ^Tr(AT^(-A^^Ae,AF, + IpJ"') + c^ 

s.t. -Tr(A^^AFj = Pr,k, 
Vk 

Y,Pr,k = Pr. (Ill) 

For any given P^ fc, then the optimization problem (lllll) can be decoupled into a collection 
of the following sub-optimization problems 

min Tr(AT,(-A»^Ae,AF, + IpJ"^) 

s.t. -Tr(AP^ AfJ = P.,fc, (112) 

''7fc 

where the constant part c^ is neglected. For any two M x M positive semi-definite Hermitian 
matrices A and B, we have Tr(AB) > ^.. Aj(A)Ajv/_j+i(B), where Aj(Z) denotes the z*'^ largest 
eigenvalue of the matrix Z ll37l . Together with the fact that elements of the diagonal matrix Atj. 
are in decreasing order, the objective function of (11121) is minimized, when {Ap A^^Ap^/rjk + 
Itv^) is a diagonal matrix with the diagonal elements in decreasing order. The objective function 
can be rewritten as 

Tr(Ai^(lAP^Ae,AF,+I^J~^) 

Ilk 

= d^(Aijd((-A«^Ae,AF, +I;vJ-^) = /(b), (113) 
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where d(Z) denotes the vector which consists of the main diagonal elements of the matrix Z. 

It follows that /(b) is a schur-concave function of b [[37l 3.H.3]. Then, based on [[TSl Theorem 
1], the optimal Ap^ has the following structure 

Apfc.opt ONk,pk-Nk 

where AF;.,opt is a A'^t x Ni^ diagonal matrix to be determined, and Nf^ = min(pfc, q^). 

Putting (11141) into the definition of Ac^^opt in (11061) . the structure of the optimal Acfc.opt is 
given by 

ACfc.opt ON,„qk-Nk 

%k-Nk,Nk ^Pk-Nk,qk-Nk 

where AG^.opt is also a Nk x Nk diagonal matrix. 
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Fig. 1. Amplify-and- forward MIMO relaying. 
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Fig. 2. MSE of received signal at the destination for different a^ when a — 0. 
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-The algorithm based on the estimated CSI only 
- The proposed robust algorithm 




E/N^ldB) 



Fig. 3. MSB of received signal at the destination for different af when ct — 0.4. 
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Fig. 4. MSE of received data at the destination for HSA, SPA and proposed robust algorithm when a = 0.4. 
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-The algorithm based on estimated CSi oniy 
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Fig. 5. MSB of received data at the destination for different a. 
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Fig. 6. BER of received data at the destination for different CTj, when a = 0.5 
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